Chameleons with Field Dependent Couplings 
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Certain scalar-tensor theories exhibit the so-called chameleon mechanism, whereby observational 
signatures of scalar fields are hidden by a combination of self-interactions and interactions with 
ambient matter. Not all scalar-tensor theories exhibit such a chameleon mechanism, which has been 
originally found in models with inverse power run-away potentials and field independent couplings 
to matter. In this paper we investigate field-theories with field- dependent couplings and a power- 
law potential for the scalar field. We show that the theory indeed is a chameleon field theory. 
We find the thin-shell solution for a spherical body and investigate the consequences for Eot-Wash 
experiments, fifth-force searches and Casimir force experiments. Requiring that the scalar-field 
evades gravitational tests, we find that the coupling is sensitive to a mass-scale which is of order of 
the Hubble scale today. 

PACS numbers: 98.80-k, 98.80.Cq, 04.50.Kd 



I. INTRODUCTION 

Modern cosmological observations strongly suggest 
that visible matter contributes only a few percent to the 
total energy budget. The rest is made of dark matter 
and dark energy. While dark matter is very well moti- 
vated within particle physics, dark energy is harder to 
explain within particle physics models. Scalar fields are 
natural candidates for dark energy, but the dark energy 
scalar field must be very light to explain the accelerated 
expansion. In addition, its coupling to matter should be 
very small. Alternatively, dark energy (and dark matter) 
might signal a breakdown of General Relativity on large 
scales. 

In the last three decades, scalar fields have played an 
important role in both cosmology and particle physics 
(see e.g. P, [H and references therein). The best mo- 
tivated particle physics candidate for a scalar field is 
the Higgs boson, part of the standard model of particle 
physics, which itself has yet to be observed. Even though 
no scalar fields have ever been observed directly yet, they 
are a general feature of high energy physics beyond the 
standard model and are often related to the presence of 
extra-dimensions. Scalar fields have been postulated as 
means to explain the early and late time acceleration of 
the Universe. However, it is almost always the case that 
such fields interact with matter: either due to a direct 
Lagrangian coupling or indirectly through a coupling to 
the Ricci scalar or as the result of quantum loop correc- 
tions. Both for inflation in the early universe and for 
dark energy, such couplings can lead to problems. In in- 
flation, for example, couplings might destroy the flatness 
of the potential needed to drive a period of inflation. If 
there are scalar fields which permeate the universe to- 
day and have non-zero couplings to matter, then they 



would induce an additional force in nature. If the scalar 
field self-interactions are negligible, then the experimen- 
tal bounds on such a field are very strong: either the 
couplings to matter are much smaller than gravity, or 
the scalar fields are very heavy, so that they mediate a 
short-ranged interaction. 

However, a certain class of theories have been pro- 
posed, in which the scalar field(s) properties depend on 
the environment: these are the class of chameleon field 
theories, proposed by Khoury and Weltman jsj, that em- 
ploys a combination of self-interaction and couplings to 
matter of the scalar-field to avoid the most restrictive of 
the current bounds. In the models that they proposed, 
which from now on will be referred to as the standard 
chameleon model (SCM), a scalar field couples to matter 
with gravitational strength, in harmony with general ex- 
pectations from string theory, whilst, at the same time, 
remaining relatively light on cosmological scales. It was 
found that local gravity constraints are (roughly) satis- 
fied as long as the mass-scale of the potential satisfies 
M < (lmm)~^. This coincides with the scale associated 
with the late time acceleration of the universe, and it is 
surprising that it should come from local experiments. 
We will, in this paper, show that this result carries over 
to other classes of chameleon models, (it will result from 
the current accuracy of the Eot-Wash and Casimir exper- 
iments), in which the coupling becomes field dependent 
and hence is environment-dependent. 

The chameleon with a constant coupling has been sub- 
ject to many studies Q-jl^l to mention some. And most 
relevant experimental bounds have been calculated for 
the two fiducial potentials introduced by Khoury and 
Weltman. There have been very few studies on the dif- 
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ferent types of coupling^. However, it would be impor- 
tant to investigate whether the chameleon mechanism is 
present in more general classes of models. Here, we will 
go one step further and generalise the chameleon mech- 
anism to an inverse power-law coupling. In doing so the 
coupling to matter becomes dynamical and will be much 
smaller on earth than in the vacuum of space. When ob- 
jects become big (in density and size) in a way defined 
later we will also have an additional suppression of the 
fifth- force by a thin-shell effect. We will show that our 
model does not yield an explicit geometrical thin-shell, as 
found by solving the field-equation in the SCM. Instead 
we find an equivalent thin-shell factor which determines 
how much the resulting fifth-force is suppressed. Adopt- 
ing the thin-shell name for our solution, we derive the 
expressions of the far-away field of thin-shelled bodies 
and show that it is independent of the parameters deter- 
mining the coupling to matter. This is the same effect as 
found in [5] for the SCM. 

This paper is divided into three main parts: in sec- 
tion II we study the behaviour of the scalar field inside 
and outside a spherical body. We find that the theory 
exhibits the chameleon mechanism and find the thin- 
shell solution. This allows us to make predictions for the 
chameleon behaviour on earth and in the solar system. 
In section HI we derive the expressions for the chameleon 
force law between different objects and ranges, which can 
be succinctly stated by introducing an effective coupling. 
In section IV we calculate the bounds on our parame- 
ters from the Eot-Wash experiment, fifth- force searches, 
post-newtonian corrections and Casimir experiments. 

We will show that the model allows for a very large 
local matter coupling, l^ph to be compatible with 
all the available data. This is entirely due to the thin- 
shell effect. We will also show how non-linear effects en- 
sure that the field value taken by the chameleon far away 
from a body with a thin-shell is independent of A, the 
parameter that describes the strength of the coupling in 
the Lagrangian. 



A. Notation and conventions 

We will always work in units of c = 1 and h = 1, 
the metric has the signature (—,+,+,+) and we will use 



the convention Mpi 



for the Planck-mass. The 



frame referring to g will be called the Einstein frame, 
and the frame referring to g the Jordan frame. When 
speaking about the chameleon mass = 14ff, </)</) we re- 
fer to the mass of oscillations about a minimum of the 
effective potential. In looking at the field inside and out- 
side a body the quantities of that body are referred to 
with a subscript c and the background with a subscript 5. 



^ See [3^ for a brief note on the power-law coupling ^(cf)) 



For example the minimum of the effective potential inside 
(outside) a body is denoted by (j)c (^b). When speaking 
about quantities such as (^5) we will sometimes sim- 
ply write 



B. The Chameleon Action 

The action governing the dynamics of a general scalar- 
tensor theory is given by 



S= dx 



V{(l>) - Cm{gtii^,1pi) 

(If 

where g is the determinant of the metric g^y^ R is the 
Ricci-scalar and ij^i are the different matter fields. The 
matter fields couple to g^y which is related to g^y via a 
conformal rescaling of the form 



(2) 



The SCM corresponds to the choice A{(j)) = e^p^ where 
/3 is a constant together with a run-away potential like 

V{(j)) = (^) ' Cosmological and local gravity ex- 



periments impose 



Mpi 



<C 1 at least since the time of 



Big Bang Nucleosynthesis (BBN) so that in most appli- 
cations of this model we can without loss of generality 
set A = 1 + -YT-- This model has been found to be 
in agreement with experiments even for /3 ^ 1 provided 
one imposes a fine-tuning in the potential no worse than a 
cosmological constant. This is different from a minimally 
coupled scalar field for which fifth-force and equivalence 
principle experiments require a coupling strength much 
smaller than unity. In this work, we will study an inverse 
power coupling 



\ogA{<P)=l3{<P) = A 



M, 



(3) 



where is a mass-scale and A a dimensionless param- 
eter. We will refer to this model as a chameleon due to 
similarities with the SCM, even though we do not know 
a priori whether this model will produce a chameleon 
thin-shell suppression effect. 



C. The Chameleon Potential 

The most important ingredient in a chameleon field 
theory is that the effective potential has a minimum 
which depends on the local matter density. The sim- 
plest type of potential, for our coupling ([3|), having this 
property is the power-law potential 



(4) 
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where M is a mass scale, n > 0, and a a dimension- 
less parameter. When n 7^ 4 we can fix a = 1 by re- 
defining the mass-scale M . This potential gives rise to 
an effective potential, defined below, of the same type 
as in the SCM. Here M can be any mass-scale, but in 
order for the chameleon to act as a dark-energy can- 
didate we need V((/)today) ^ A together with an equa- 
tion of state uo ^ —1. It is therefore convenient to set 
M = M]jE = A4 and have the cosmological constant as 
part of the potential. In this case we can think of the 
potential as a Taylor expansion of a more complicated 
potential such as V = M"^ exp((^^/M^), for <C M. 



D. The Field equation 

Variation of the action ([1]) with respect to (/> yields the 
field-equation 



2 dCr 



(5) 



where the sum is over the different matter-species and we 
have allowed for different couplings to different species. 
Assuming that the matter fields tjj'^ do not interact with 
each other, each energy- momentum tensor (suppressing 
the (i) for now) 



2 dCr^ 



is conserved in the Jordan-frame l39| 



(6) 



is the Einstein-frame density. With this definition, the 
equation of motion, in the Einstein- frame, reads 



(l-3u;i) 



(12) 



and we see that the field equation for (j) can be written 
□(/) = 14ff,(/)7 where the effective potential is given by 

FeffW = + ^Mf'^^^^W (13) 



To simplify things we will assume that all the different 
matter species couple to (j) with the same A{(j)) and we 
will only consider non-relativistic matter where uJi ^ 
With these assumptions, the effective potential reduces 
to 



(14) 



Note that since the matter fields couple to ^, the 
geodesies of a test-particle will be the geodesies of this 
metric, and p is the physical density. We do not need 
to be too careful about this since, as we will show, in all 
practical applications we will have A{(j)) ^ 1 and the two 
densities are essentially the same. 



E. Minimum of the effective potential 

The minimum of the effective potential is determined 
by the equation ]4ff,0 = which gives 



(7) 



where V is the Levi-Civita connection corresponding to 
the metric g. In the perfect fluid approximation where 
each matter species behaves as a perfect isentropic fluid 
with equation of state p = ujip we have 



(8) 



Going to the Einstein frame we choose, without loss 
of generality, a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) background metric. The energy density p in the 
Einstein-frame is the one that obeys the usual continu- 
ity equation p oc a~^^^~^^^\ Computing the Christoffel- 
symbol 



_ /AM/3V+^ / kp 



The chameleon mass at the minimum is given by 



, N fe(n+fe)p / AM^\ 



k+2 



(15) 



(16) 



where we have used that in contrast to the standard 
chameleon where m| = F,^^ we now have to take in ac- 
count the contribution from the term p. But we can 
ignore the term pP'^^c/) which is valid as long as /3{(j)) <C 1. 
From (p!6|) we see that the field can only be a chameleon 
for n > 2. 



and using ([7j) we find 



Thus 



(10) 



(11) 



F. An equivalent formulation 



We redefine the field by introducing x = ~^^v>i' The 
coupling yields 



(17) 
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which becomes that of the SCM for k 
law potential (|4]) becomes 



vix) 



4 f 



\ X 



1. Our power- 



(18) 



where M« = When = ^ we have = 

M and the potential is identical to the Rattra-Peebles 
potential often used in the SCM. With this choice for 
M/3, the full action can be written 



S = / dx^^f^ 



(19) 



and we see that (for /c = 1) it is only the kinetic terms 
that distinguish our model from the SCM. The fine- 
tuning in the coupling sector is removed and we are left 
with only one fine-tuned mass-scale in the action. 
The field equation is given by 



□x-|(V,x)' = (^)*T4ff,x^ 



(20) 



which is significantly more complicated to work with than 
(|5]) so we will use the original formulation. 



G. The Coupling Scale 



II. SPHERICAL SOLUTIONS TO THE FIELD 
EQUATION 

The field equation in a static spherical symmetric met- 
ric with weak gravity reads 



r 



(25) 



where we have assumed I3{(j)) <C 1. We study solutions 
inside and outside a spherical body of constant density 
Pc (e.g. the earth) in a background of a very low density 
Pb < Pc, 



pc for r < R 
Pb for r > R 

and impose the boundary conditions 



-{ 



(26) 



d4 
dr 



dr 



r=0 



= 
= 



(27) 



The first condition follows from the spherical symmetry 
around r = and the second one implies that the field 
converges to the minimum of the effective potential, 
in the far-away background. If rricR <C 1, the chameleon 
acts approximately as a linear scalar field whereas in the 
case rricR ^ 1 the full non-linearity of the field equation 
comes into play. 



In the background today, taking A = 1, we have 
(t)b /M\^/Mde\^ 



(21) 



For the model to be in agreement with experiments we 
must require /3{(j)i)) ^ 1- This constrains 



Mf3 <C Mde 



M 



(22) 



showing the need to fine-tune Mp. We fix by the 
requirement that the equivalent action (|T9l) is of the same 
form as the SCM when M = Mde ~ {lmm)~^. This 
fixes 



M« = 



Mpi 



(23) 



This choice also ensures that the coupling |Mpi of a 
pc ^ Ig/cm^ body is of order 1 when M = Mde- The 
term \Mp\ plays the same role in this model as (3 does 
in the SCM, but here this factor is field-dependent. In 
the remaining of this article we take = Hq so that 



A- 



(24) 



is our coupling. 



A. Case 1: the Thick-shell rricR <C 1 

In this case the initial value satisfies (j){0) = <^ 
(j)c and the approximation ]4ff,(/) ^ P^cf) Pc is valid inside 
the body. Since this driving force is relatively small, we 
approximate /3,0Pc ^ /^,0zPc- Solving the field equation is 
now straightforward and the solution reads 



6 



for < r < 



(28) 



where we have used absolute values since < 0. This 
solution corresponds to the thick-shell solution in the 
SCM, which is not surprising since the non-linearities in 
the field equation are negligible. Outside the body we 
assume that the linear approximation is valid leading to 
a Yukawa profile (/) = -\- — Matching the two 
solutions ai r = R leads to 



with determined through 



(29) 



(30) 
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Defining mf = pcP^cfycf) {4^i)i tiie cilameleon mass at ttie 
centre of tiie object, tiiis last expression can be rewritten 
as 



(31) 



and the chameleon takes a value at the centre of the body 
corresponding to a mass ~ This also shows that 
the approximation used inside the body is valid since the 
field undergoes a "^^^^"^^^^^ < (9(1) change. The initial 
value, can be rewritten in a more compact fashion for 
> as 



(32) 



If nicR is really small we have ~ and the field inside 
the body is just a small perturbation in the background. 
To summarise, the solution is 



= ^6 + 
= ^6 + 



per 
6 

1/3, <^J Mie' 



47r 



for <r < R 
foT R < r 



(33) 



The far-away field is proportional to the coupling 
evaluated inside the body (or equivalently at the surface), 
just like for the SCM. Let us mention that two bodies 
with rricR <C 1 attract each other with an attractive force 
of magnitude 



(2) 



GMiM2{l-\-mbr)e- 



(34) 



The relative strength to gravity can be read off as 
2/3,^(1) /3,^(2) which is maximal for bodies where ^ 

^5. If a body increases in size the strength of the fifth- 
force decreases. In contrast with the SCM, this sup- 
pression appears even for bodies without a thin-shell 
{rricR <C 1). See Fig. [T]for a plot of a rricR <C 1 profile 
compared to the analytical approximation found above. 

Note that outside the body we have assumed that the 
Yukawa profile is a good approximation. When 
05 we have ^(i^) ^ ^5 meaning that the approximation 
Vei^^cf) ^ ^li^j) — 4^b) is not valid right outside r = R. In 
these cases, the driving term y,^ can be neglected relative 
to the friction term, leading to the same 1 /r-profile. This 
approximation is valid up to the region where m^r ^ 1 or 
equivalently (j) ^ which leads to the Yukawa solution 
used above. As rribR <C 1, we can add the exponential 
factor to the solution outside r = R. The numerical 
results show that the analytical solutions found above 
match the actual solutions to a good level of accuracy. 



B. Case 2: the thin-shell rricR > 1 

In this case the field starts out very close to the min- 
imum inside the body: ~ (j)c- The field will re- 
main close to the minimum throughout the body making 



1. 
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FIG. 1. Numerical field profile for {n,k) — (10,1) and 
rricR = 10 ~^ together with the analytical approximation 
(dashed line). The analytical approximation is seen to be 
a very good match to the actual solution. 



the linear approximation V^ff,^ = m^(0 — (j)c) valid in 
< r < R. This case is similar to the thin-shell solution 
in the SCM, with the exception that we do not have this 
explicit thin-shell. Nevertheless we will adopt the thin- 
shell name for our solutions. Right outside a thin-shelled 
body the approximation 14ff,0 ~ y,^ is valid and we must 
solve 



r 



naM^ i- 
M 



for < r < i^* 



(35) 



where R* is the point where the coupling term, pbP^^, 
becomes relevant again. When n = 1 or n = 2 we 
can solve (|35|) as it stands. In these cases the field is 
not a chameleon since the chameleon mass increases as 
the density decreases. In the general case we will need 
certain approximations to find a solution. 

In Fig. [2] we plot a thin-shelled solution for the 
earth in the cosmological background (density equal to 
the average cosmological density) for n = 10. 




FIG. 2. The Thin-shell profile for the earth when (n, k) = 
(10,1) and mcR= 10^ 

Inside the body the field is very close to the minimum 



6 



and remains there throughout most of the body, except 
near the surface where the field undergoes a small change. 
Linearizing the effective potential around 0c- Kiff,^ = 
m^(0 — 0c) 7 we can find the solution close to r = that 
matches the initial condition 



= 00 (l-r^i^) 



in < r < i? 



4>c-4>{0) 



< 1. 



(36) 



The solution is valid as long as the linear term in the 
Taylor expansion of 14ff,0 dominates over the higher order 
terms, which gives the condition 



< 



< 



|n-fe-3| 
((/c+lK/c+2)) 



1/2 



for n — /c — 3 7^ 
for n — /c — 3 = 0. 



(37) 



The largest value of |0 — 0c| inside the body occurs at 
r = R and we will later check that this value satisfies the 
condition above. Defining S = ^EBMii^^ haxe that 
the field value and derivative at r = ii satisfy 



(1 - S)^c 

-6mc(l)c 



(38) 



Outside the body the potential is very steep, so the field 
drops very quickly and the friction term ^0 is, initially, 
negligible compared to the driving force F,^, implying 
that 



To simplify the analysis we define ij) 



= ± 



(39) 
and 



— =' so that we can write the equations in a dimension- 



dx 

less form as 



which has the solution 

i, = 
a 



[l+a(r/i?-l)]— 2 



y'2n(n+/c) 

Matching to the solution for r < using 

_ 2 
(l-(5)^ ~ n(n + /c) 



(40) 

(41) 

we find 
(42) 



which determine^ S. Numerically we find S ^ 0.086 
when (n, k) = (10, 1) and S ^ 0.200 when (n, k) = (4, 1) 



^ The reader may question this result, since when 6 is determined 
the full solution in < r < oo is known (at least from a nu- 
merical point of view), but it is derived without considering the 
behaviour at large r yet. What this result really states is that 
the solution that converges to in the background will have to 
correspond to a particular initial value of order 6 at r = R. This 
result is confirmed by the numerics. 



for all rricR > 10 in very good agreement with the formula 
above. If we now go back and put this value for 6 into 
(|37|) we see that the linear approximation is valid for all 
reasonable values of {n^k). As an example, take n = 4, 



upon using ([32]) we find 8'^ 
the condition 



1 



2(/c+4) 



which in (|37|) gives 



1 < 



8(fe + 4) 
{k-lf 



k < 12. 



(43) 



As the field rolls down along the potential, it reaches a 
point r = R^ where the driving force satisfied 



rivmg 



< 1 



(44) 



and from here on the dynamics of t/j are determined by 
the friction term which we have neglected until now. The 
field equation reads 



+ ~ for i?* < r 

X 

with the solution 

il; ^ -\- for R"" < r 



(45) 



(46) 



for some A. This solution is valid until we reach the re- 
gion where the driving force has to be taken into account 
again. This is the case when m^r ~ 1 or equivalently 
i/j ^ i/ji) and alters the solution by adding a Yukawa ex- 
ponential e~^^^ to the 1/r term. Again since m^i?* < 1 
we can incorporate this by adding this term to (|46|) as 



^^*g-m,(r-i?*) 



for i?* < r < OO 



(47) 



The matching of (|4T]) and (j471) at r = i?*, defining A = 
^ implies the identifications 



R 



A 



(l+aA)'^-2 

(1+A)a 2 



(48) 



(l+aA)^ 



1+aA n-2 



When mi)R < 1 we have ip{R*) 
leads to 



A -\- t/ji) ^ ijjh which 



A 

AR* 

B 



BR 



BR 



{nicRy 

( n(n+fe) \^ / n-2\ 
- \ 2 ) \n-A) 



{rribR) ' 



(49) 



where we have used aA 1 in order to simplify the 
solutions. This derivation does not apply for n = 4. A 



^ In the case where m^i? ^ 1 the field will have settled at the 
minimum before this happens. Since this corresponds to a very 
short ranged force we will not consider it here. 
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similar derivation shows that (|49|) is vahd for n = 4 when 
one takes the hmit n ^ 4 in the expression for B. Let us 
summarise the solutions we have found: 



for r < R 

for R<r < R* 



(l+a(r/i?-l))— 2 

06+ ^e-^"--"'' for R*<r. 



(50) 



Defining the effective coupling in the thin-shell case via 



we have that 



47rM 



for i?* < 



pi 



n-2 
n - 4 



(2a)- 



(51) 



(52) 



which is independent of the parameters defining the cou- 
pling. Thus, the exterior profile of a thin-shelled body 
depends only on the radii R^ the mass Mi and the po- 
tential parameters. In the SCM, using the same potential 
as we do, the same effect was found in [5]. Our expres- 
sion for the effective coupling ([52|) agrees with [B', Eq. 
24]. This similarity should not come as a surprise since 
a thin-shell solution is associated with dominating self- 
interactions. 

To justify the thin-shell name for our solutions, we 
define a thin-shell factor via 



/3eff= |/3,^jMpi 



and we find 



B 



\P.^JpcR\mcR)^ {rricR) 



2(r^-l) 



(53) 



(54) 



This factor determines how much of the mass of the body 
contributes to the fifth- force. As rUcR ^ 1, we have 
<C 1 and thus /3eff <C l/^,^^ l^pi- If we extend this 
definition and set /3eff = |Mpi when rricR <^ I then 
(|5T]) is valid for all bodies. See Fig. [3] for a plot of the 
effective coupling as a function of the radius of the body. 
We note that for the special case when n = 4 the far-away 
field can be written as 



rrib 



for R* <r 



(55) 



which is completely independent of the parameters pc and 
R describing the body, and depends only on the back- 
ground. Likewise /3eff only depends on the mass of the 
body. We are interested in the cases where the field has a 
large range in the solar-system together with thin-shelled 

n— 2 n—2 

planets m^R » 1. From ^ = {^) ' = (^) 

we see that having a small k and a large n gives the 







— 


^i = 0l. / 








-10 




-u 




■ / \ 





FIG. 3. The effective coupling for a spherical body with con- 
stant density. When the body is very small the field inside 
the body is the same as the background, ^b, leading to a big 
coupling. Then as the radius gets bigger the field inside the 
body starts moving away from the background and the cou- 
pling decreases. Finally when we reach rricR > 1, the field 
inside the body settles at (/)c and develops a thin-shell such 
that the coupling starts to decrease like 1/R^. 



largest ratio ^ ~ ( — I • For the case of the earth 

pc ~ Ig/cm? in a background of the average solar sys- 
tem density ph ~ 10~^^^/cm^ we find ^ ~ 10"^^. It is 
therefore possible for the field to have a range as large as 



10 



15^ 



l{)^Au and at the same time to have a 



thin-shelled earth rricR ^ 1. This is the same as found 
in the SCM 

In the SCM, the coupling is easily identified as the pa- 
rameter that multiplies (j) in the matter-Lagrangian. Here 
we have a coupling that varies from place to place and is 
in general given by /3eff defined above. A test particle in a 
region where (j) ^ (j)o will experience a coupling |/3,0q \Mp\. 
The coupling becomes smaller in a high-density environ- 
ment and the largest value is achieved in the cosmolog- 
ical background. One can say that the chameleon effect 
in this model is twofold: first the coupling decreases as 
the environment gets denser and secondly for really large 
objects only a thin-shell near the surface contributes to 
the fifth force. 

In Fig. m we see the numerical thin-shell profile together 
with the analytical approximation for (n, k) = (10, 1) and 
rricR = 10^. To calculate the profile for such highly thin- 
shelled objects {rricR > 100) it is not possible to start 
the numerical simulation at r = since the initial value 
is too close to (j)c- Upon using the relation (f38|) and (f42|) 
we are able to start the simulation at r = allowing us 
to produce the field profiles shown here. 

Note that for large rricR we have a large gradient at 
r = R that may cause problems in laboratory experi- 
ments using a very small separation between objects (like 
Casimir, Eot-Wash etc.). 
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FIG. 4. The thin-shell profile for the earth when (n, k) — 
(10, 1) and TTicR — 10^ together with the analytical approxi- 
mation (dashed line). The horizontal line shows the min- 
imum in the background. The error between the numerical 
solution and the analytical approximation is less than 10% in 
the whole range. 



III. THE CHAMELEON FORCE 



The geodesic equation in the Jordan frame reads 



(56) 



Using this can be rewritten in terms of the Einstein 
frame connection T and 6 as 



(57) 



In the non-relativistic limit the last term can be neglected 
and the chameleon force on a test particle is given by 



(58) 



This is attractive since both and are negative 
outside a spherical object as shown in section IL 



A. Chameleonic Force between two parallel plates 

We consider the force between two identical parallel 
plates of radius Rp whose surfaces are separated by a 
distance d <^ Rp and the system is in a laboratory vac- 
uum (5|,[i8|. In practice the 'vacuum' will have a non zero 
pressure corresponding to a very low, but non-zero den- 
sity pb. Because the plates are very close to each other 
we can treat the plates as infinite flat slabs and take plate 
1 to occupy the region x < —d/2 and plate 2 to occupy 
the region x > d/2. 

Let us first consider the case when the scalar field is 
linear with a constant mass m and coupling A. Since lin- 
earity means the superposition principle holds, we only 
need to consider the field emanating from plate 1 in or- 
der to calculate the force between the plates. The field- 
equation for a linear scalar field reads 



m^(p + — — 
Mpi 



(59) 



where = because of the symmetry in the setup. 
The field equation has the solution 



Ape 



for X < 



5e-^^ for X > 



(60) 
(61) 



Matcfiing tfie two solutions at x = gives us B = —C = 
^^"^ ~ Tfie force on tfie second plate due to tfie first 



2M„ 



is then given by 
A 



Mpi 

SttA' 



f 



\/(j)dx 



Ape 



(62) 
(63) 



where A is the surface area of the plates. Returning to 
our model, we use a subscript s when talking about the 
quantities defined at the surface of the plates, subscript h 
in the background and subscript c inside the plates. For 
example the field- value at the surface of the plates is re- 
ferred to as (j)s^ y{(l)c) = Vc and so on. Also a subscript 
is used to refer to the quantities where (j) = between 
the plates. Because of the symmetry this is at the point 
X = 0. Finally we assume that the chameleon mass sat- 
isfies rricRp ^ 1 so that the true non-linear nature of the 
chameleon comes into play. With the conditions stated 
at the beginning, we have that (j) obeys 



d^(^ 



between the plates, and 

^ _ 
dx'^ 



(64) 



(65) 



inside either plate. Integrating the equations above yields 



= 2{V{4>) 
'f = 2(y((/.) 



Pb(/3(0) - /3o)) for -dl2<x< d/2 
■ Peim - Pc)) forx2>dV4. 

(66) 



Where we have used that deep inside the plates (/)(±(X)) ^ 
and ^^^^ = 0. Matching at x = ±d/2 we find that 



the coupling at the surface is given by 



Vc-Vo+ PcPc - PbPo 
Pc- Pb 



(67) 



If the second plate were removed (j)o = (pb^ the coupling 
at the surface /3so would be given by (|67|) . with (j)o ^ (j)}). 
The perturbation, S/3s = Ps — /^so, in P{<Ps) due to the 
presence of the second plate is therefore 



Vb-Vo + Pb{lib-M 



(68) 



where we have used pc ^ ph- Since rricRp ^ 1 the per- 
turbation deep inside the plates are suppressed exponen- 
tially. Using (1) we find that the attractive force on one 
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plate due to the presence of the other one is given by 



\3t 

A 



Id/2 dx 



is exponentiahy suppressed. We put 0o = ^6(1 + ^) 
where we assume (5 <C 1. This allows us to approximate 

n-\-k 

^ 1 — (n + and since (j)c ^ </>6 we can take 
^5 ^ oo and the integral ([72]) can be written 



which, using ([67|) . gives 

^ = ^0 - n + P6(/3o - M = K.ff((/>o) - K.ff((/>6) (TO) 

We have to calculate the field value (pQ halfway between 
the plates. This is done by integrating (|66|) over the 
region —d/2 < x < 0, using that ^ < when taking the 
square root. This gives the equation for 



(71) 



dz 



rribd 



y'^" - 1 + f (1 - (n + k)5) - 1) \/2n(n + k) 

(76) 

In the limit (5^0 the left hand side diverges. Upon using 
a power series expansion of the integrand near z = 1 

^"-l + f(|^) (^-'=-l)«n(n + fc)x (77) 
[(5(^-l) + i(l-(5(fc+l))(^-l)2 



00 VV{4>)-Vo + Pb{l3{4>)-Po) V2 



(n-fc-3)-(fc+l)(fc+2)(5 



(^-1)3 + ... 



(78) 



This is a general expression, and can be used for any 
coupling and potential. Specialising to our case where 

m = (a^)' and V{^) = cTM^(^^y^ we change 

variables to z = (j)/(j)o and define Zs = 0s /^o giving 



we see that the second term is the divergent part when 
^ = 0. This term dominates in the region 1 -\- 2S < z < 

and for < J <C 1 provides the dominating con- 



I n-k-3 

tribution to the integral. We can therefore approximate 
the integral by 



i: 



dz 



(72) 



Here we can have several cases. 

1. Case 1: 00 ~ 0c 



This case corresponds to very small separations rrihd <C 
1. We set 00 = 0c(l — and rewrite the right hand side 



I n-k-3 I 



dz 



/1+25 

This gives 



1 rn^d 

and shows that the chameleon force 



A 



eff,< 



(^0 - 4>b)' 



(79) 



(80) 



(81) 



of Cl as 



-y/2n(n+fe) 



. The integral can now be evaluated jg indeed exponentially suppressed by the factor rribd > 



1. 



r 



resulting in 



dz 2^/^ 

^n{z - 1) 



{rricdf 



(73) 



(74) 



Case 3: 



> 00 > 06 



In this last case we can neglect the third term in the 
square root of ([72|) and also take ^ oo. This enables 
us to evaluate the integral analytically 



This case only applies when the separation d is much 
smaller than the thickness of the plates t, since we have 
assumed rrict ^ 1. The chameleon force becomes 



dz _r(i)r(i 



|r(-^)l 



(82) 



A 



n{mcd)'^ 
8{n + k) 



(75) 



2. Case 2: 0o ^ 06 

This corresponds to the case when the field drops all 
the way down to the minimum in between the bodies. 
Since this case corresponds to m^i? > mi^d > 1 the force 



The F-function satisfies r(e) ^ — 7^; for e <C 1 with 
7£; ^ 0.577 being the Euler- Gamma constant. This gives 

= ^^tw ^ - i^^g^ ^- (S3) 

|r(-^)l ri 

We can now find an explicit expression for 0o 



(84) 
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Beha™ur of Chameleonic Pressure 



Constant 
i-^--^ j- Force 
f> Behavior 

< 




Separation of plates: d 

FIG. 5. General behaviour of the chameleon pressure as 
a function of the plate separation d. 



and the chameleon force 



2 Sfi 
a Md 



(85) 



We see that the force follows a power-law where the drop- 
off is faster than l/d^, but slower than l/d^ when n > 4. 
The Casimir- force falls off as 1/d'^ making Casimir ex- 
periments (with large plate-separations) a powerful way 
of constraining the chameleon. See Fig. [5] for a plot of 
the chameleonic force (or more accurately, the pressure 
F(f)/A) as a function of the distance between the plates. 



Chameleon Force between two spherical 
thin-shelled bodies 



far field of body one evaluated at body two 



47rM 



pi 



(87) 



Because of the big mass of the chameleon inside body 
two, the perturbation created by body one is attenuated, 
and, as in the SCM, only a thin-shell close to the surface 
contributes to the force. We model this by setting 

\F\ 9R B M (^^\ GMMl+m,r)e-"^^- 



R J. 



(88) 



where (•^)2 models the effect of this thin-shell. 

Likewise the force on body one due to body two is 
given by the same expression with 1 — > 2. Up to a 0{\) 
factor we have 



R 



(89) 



which equals |^ ^""^j^ ^ up to a 0{1) factor. The force 
between two thin-shelled objects is then given by 

|i^0| = ^PefflPeff2 ^ (90) 

where we have chosen an appropriate 0{1) factor. In the 
thick-shell case {rricR <C 1) the whole body contributes 
to the forccH 

iF,i=2<«„)(ftg'M.,) °"-"^'":'""-'-'"" 

(91) 



We consider the force between two bodies, with thin- 
shells, that are separated by a distance r i?i,i?2- 
Given than r ^ Ri^R2 we can consider the monopole 
moment of the field emanating from the two bodies only. 

We denote by {(j)2) the field outside body one (two) 
when body two (one) is absent. To a good accuracy we 
have ^ (pel the minimum inside body 1. In between 
the bodies, we can superimpose the far-away fields from 
the two bodies. As the distance is large, the perturbation 
5(j)i in the field inside the body two due to the presence 
of body one satisfies 5(j)i <^ 02 • The combined field close 
to the surface of body two is approximately given by 

(j)2 + ^01. 

Using the geodesic equation, dF(f) = — \/(j)dm, we 
have that the total force on body two due to body one is 



IV. BOUNDS ON THE PARAMETERS 

We will constrain the parameters A and M (or a) by 
looking at the consequences our model has on local grav- 
ity experiments. The experiments considered here re- 
stricts the value of the chameleon coupling in different 
regions. The Eot-Wash experiment (and other fifth-force 
searches) are usually the best way to obtain good bounds 
when 1/3,0^ |Mpi ~ 1. Casimir type experiments are often 
the best way to bound the highly coupled, |Mpi 1, 
region. Finally the PPN and BBN bounds constrain 
the extremely high coupled region which are invisible to 
the Casimir type experiments due to the extremely short 
range of the chameleon. 



1^01^/3,0,, / \/6(i)idm (86) 

^body two 



Next we have that the perturbation 5(j)i is given by the 



The field-equation is quasi-linear and the superposition principle 
holds. 
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A. PPN bounds 

For experiments using the deflection of light by large 
bodies, the only Post-Newtonian Parameter (PPN) at 
play is the Eddington-parameter 7. The Eddington- 
parameter is defined in the Jordan- frame by gij = (1 — 
2'y^)Sij when ^00 = — 1 — 2^ [41]. Transforming to the 
Einstein-frame we get the following estimate for 7 



PPN bounds oiiiif=4. k^I liieories 



7 



1 



(92) 



The back reaction on the gravitational potential from 
the chameleon is in most interesting cases negligible, and 
since ^{(j)) <C 1 the Jordan-frame and Einstein-frame po- 
tential are the same. The best bounds on this param- 
eter comes from the Cassini-experiment (42| and reads 
I7 — 1| < 2.3 • 10~^. The gravitational potential for the 
sun is ^sun = 10~^ and the field near the surface of the 
sun satisfies ^ 0c giving us the bounds shown in Fig. [6l 
This experiment only restricts the parameters in which 
1/3,^, |Mpi » 1. 



B. BBN bounds 

Since our chameleon couples to matter via the con- 
formal transformation ([2]), the masses of the standard 
model particles have a (/)-dependence of the form m = 
moexp/3(0). Bounds on particle masses restrict a varia- 
tion of this type to be below the 10% level since Big-Bang 
Nucleosynthesis (BBN) Q. Since in our model (/> < 0, 
I3{(j)) is an increasing function of time so we must require 



/3((/>today) < 0.1 and /3((/)bbn) < 0.1 



(93) 



The last condition is satisfied as long as the chameleon 
has settled at the minimum before the time of BBN. The 
condition today translates into the bound 



\MdeJ 
A < lO^V-^ forn = 4 



(95) 



The k dependence is weak, and we have that this bound 
is satisfied as long as the PPN bound above is satisfied. 



C. Eot-Wash bounds 

The University of Washington Eot-Wash experiment 
(4Q| is designed to search for deviations from the 1/r^ 
drop-off of Newton's law. The experiment uses a rotat- 
ing torsion balance to measure the torque on a pendu- 
lum. The torque on the pendulum is induced by an 
attractor which rotates with a frequency uj. The at- 
tractor has 42 equally spaced holes, or missing masses. 
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FIG. 6. PPN constraints on chameleon theories coming from 
experimental bounds on the Eddington-parameter in light- 
deflection experiments. The shaded area shows the regions of 
parameter space that are allowed by the current data. The 
solid horizontal black lines indicate the cases where M and 
a take 'natural values'. The solid vertical lines show when 
M/3 = Hq. The dashed black line indicates when |Mpi = 
1 for pc = 0(lg/cm^). The amount of allowed parameter 
space increases with n. 
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bored into it. As a result, any torque on the pendu- 
lum, which is produced by the attractor, will have a 
characteristic frequency which is some integer multiple 
of 21uj. This characteristic frequency allows any torque 
due to background forces to be identified. The attrac- 
tor is manufactured so that, if gravity drops off as 1/r^, 
the torque on the pendulum vanishes. The experiment 
has been running with different separations between the 
pendulum and the attractor. The experiment has been 
running for typically separations d = bbfim. Both the 
attractor and the pendulum are made out of molybde- 
num with a density of about pc = lOg/cm^ and are 
t = 0.997mm thick. Electrostatic forces are shielded 
by placing a (ishieid = lOfim thick, uniform BeCu sheet 
between the attractor and the pendulum. The density 
of this sheet is Pshieid = SAg/cm^. As discussed in ^] 
the role played by this sheet is crucial when testing for 
chameleon fields in the strong coupling regime. If the 
coupling is strong enough, the sheet will itself develop a 
thin-shell. When this occurs the effect of the sheet is not 
only to shield electrostatic forces, but also to block any 
chameleon force originating from the attractor. Follow- 
ing the analogy of our model with the SCM this effect 
is given by an extra suppression of e"^^^^^^^^^^^^^^. And, 
in effect, this will make a larger part of the parameter 
space allowed in the strongly coupled case. It will not 
affect the experiment when |Mpi ^ 1. The force per 
unit area between the attractor and the pendulum plates 
due to a scalar field with matter coupling A and constant 
mass m, where 1/m <C 0.997mm is given by (|62|) 



A 



= a- 



Gpl 



2m? 



(96) 



where a = SttA^ and d is the separation of the two plates. 
The strongest bound on a coming from the Eot-Wash 
experiment is a < 2.5 • 10~^ for 1/m = 0.4 — O.Smm 
which constrains A < 10~^. 

When the pendulum and the attractor have thin- 
shells the force is given by the expressions derived in 
section III. The vacuum used in these experiments has 
a pressure of p = 10~^ Torr which means that the 
chameleon mass in the background, m^, is non-zero 
and for the largest couplings we will have a e~^^^ 
suppression. Hence the experiment cannot detect a very 
strongly coupled chameleon. The BeCu sheet produces 
a force on the pendulum. As the sheet is uniform, this 
resulting force leads to no detectable torque. If neither 
the pendulum nor the attractor have thin-shells then 
we must have mijd <C 1 and the chameleon force is 
simply 2/3,^. times the gravitational one. Since this 
force drops off as 1/r^, it will be undetectable in this 
experiment. In this case, however, the model-parameters 
are constrained by other experiments such as those that 
look for Yukawa forces with larger ranges as discussed 
below. 

Even though we have formulae for the force, we have 
used numerics to calculate the bounds. This gives more 



accuracy in the regions where our approximate formulae 
do not apply. The torque generated by the rotation of 
the plates can be shown to be given by [18] 



r0 



-^shieldC^shi€ 



dx 



(97) 



where ax = ^ is a constant that depends on the experi- 
mental setup and the exponential models the effect of the 
electrostatic shield. For the 2006 Eot-Wash experiment 
ttT = 3 • 10~^m^. The bounds derived from the exper- 
iment can also be expressed in terms of this torque as 
T(i){d = bbjim) < 0.87 • 10~^^Nm, which we have used to 
compute the bounds numerically. We have also compared 
the numerical results and the analytical expression in the 
regions where they both apply. Our results are shown in 
Fig. In these plots the shaded region is allowed by the 
current bounds. 

When n = 4, we can see that a natural value of a is 
ruled out for A = 1. As n becomes larger than 8, the case 
A ~ 1 becomes allowed for M = Mde- There exists, for 
each n, a large region of the parameter space which is 
allowed by the experiment and in which Mp\ ^ 1. 

The area of allowed parameter space grows with in- 
creasing n. Indeed when the potential is steeper, the 
mass of the chameleon increases, and the thin-shell effect 
is present for a larger part of the parameter space. 

The setup and the behaviour of a chameleon in the 
experiment is more thoroughly explained in [l8|. 



D. Fifth-force searches 

In the Irvine-experiment (43| the inverse-square dis- 
tance dependence of the Newtonian gravitational force 
law was tested. One experiment used a torsion balance 
consisting of a 60-cm-long copper bar suspended at its 
midpoint by a tungsten wire, to compare the torque pro- 
duced by copper masses 105 cm from the balance axis 
with the torque produced by a copper mass 5 cm from 
the side of the balance bar, near its end. The produced 
torques due to the masses at 105 cm and 5 cm have been 
measured. Letting i^Measured be the measured ratio of 
the two torques and i^Newton the Newtonian prediction 
it was found that 



-^Measured 



(1.2 ±7) • 10" 



(98) 



If the walls of the vacuum chamber do not have thin- 
shells, the field inside the chamber, as discussed below 
(jSl]), settles at a value where mchamber ~ ^^limber where 
^chamber IS the sizc of the chamber. The experiment here 
bounds 



2/3,^(1) Ml < 10-^ 



(99) 



with determined by (|3Q|) . The vacuum chamber used 
was held at a pressure p = 3 • 10~^torr which corresponds 
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FIG. 7. Constraints on chameleon theories coming from Eot- 
Wash bounds on deviations from Newton's law. The shaded 
area shows the regions of parameter space that are allowed by 
the current data. The solid horizontal black lines indicate the 
cases where M and a take 'natural values'. The solid vertical 
lines show when — Hq. The dashed black line indicates 
when 1/3,0^ |Mpi = 1 for pc = 0(lg/cm^). The amount of 
allowed parameter space increases with n. 



to a background density 4.6- 10~^^^/cm^ (at T = 300K). 
When the walls of the chamber (and therefore the test- 
masses) have thin-shells the chameleon sits at the min- 
imum of its effective potential inside the chamber. The 
chameleon mass mchamber is typically much larger than 
the inverse size of the chamber and the bound becomes 



2/3efn/3eff2 < 10- 



(100) 



where /3eff is the thin-shell effective coupling given by 
([52|). For the highly coupled cases mchamber Rchamber > 1 



there is an extra e"^^^^"^'^®^^ suppression of the torque 
where d is the separation of the test-masses. This exper- 
iment provides the best bounds for the chameleon in the 
linear regime since the more accurate Eot-Wash experi- 
ment is, by design, unable to detect the linear chameleon 
{F(jy (X 1/r^). See Fig. [8] for the resulting bounds. 



E. Casimir bounds 

Casimir force experiments provide an excellent way 
of bounding chameleon field parameters when the scalar 
field is strongly coupled to matter. Casimir force exper- 
iments measure the force per unit area between two test 
masses separated by a distance d. It is generally the case 
that d is small compared to the curvature of the surface 
of the two bodies and so the test masses can be modeled, 
to a good approximation, as flat plates and the results 
derived in section III apply. The Casimir force between 
two parallel plates is: 



li^Ca 



240^4 



(101) 



Even though the most accurate measurements of the 
Casimir force have been made using one sphere and one 
slab as the test bodies, this setup has a more complicated 
geometry and will not be discussed in this paper. We will 
focus on the experiments which use two flat slabs as test 
bodies. 

In all cases, apart from n = 4 and rricd ^ 1, the 
chameleon force per area grows more slowly than d'^ as 
d ^ 0. When n = 4 and rricd ^ 1, m^jd <C 1 we 
have F(f)/A ex d"^. It follows that the larger the sep- 
aration, the better Casimir force searches constrain 
chameleon theories. Additionally, these tests provide the 
best bounds when the test masses do have thin-shells as 
this results in a strongly d dependent chameleon force. 

Note that if the background chameleon mass is large 
enough m^d ^ 1 then i^^/A is suppressed by a factor 
g-mbd This shows that the experiments cannot detect 
the strongest coupled chameleons. For these extreme 
cases the post-newtonian corrections (and BBN bounds) 
constrain these theories. See (i^l for a detailed analysis 
of the Casimir force in the SCM. 

To date, the most accurate measurements of the 
Casimir force over separations d = 0.16 — 1.2/im have 
been made by Decca et al. in a series of three experi- 
ments taking place between 2003 and 2007 (i5|-[i3- We 
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FIG. 8. Constraints on chameleon theories coming from ex- 
perimental fifth- force searches (the Irvine-experiment). The 
shaded area shows the regions of parameter space that are 
allowed by the current data. The solid horizontal black lines 
indicate the cases where M and a take 'natural values'. The 
solid vertical lines show when — Hq. The dashed black 
line indicates when 1/3,0^ \Mpi — 1 for pc — 0(lg / cm^). 



define P = ^ to be the total measured pressure between 
two parallel plates. Using their most recent experiment, 
described in Ref. ji^l, Decca et al. found the follow- 
ing 95% confidence intervals on AP = P — Pcasimir* 
at d = 162nm, |AP| < 21.2mPa, at d = 400nm, 
|AP| < 0.69mPa and at d = 746nm, |AP| < 0.35mPa. 
The resulting bounds are shown in Fig. [9l The area of 
allowed parameter space grows with n and k. 



F. Combined bounds 

The chameleon theories considered in this work have 
a four-dimensional parameter space, spanned either by 
M and A (n > 4), or by a and A (n = 4). We combine 
the constraints found in sections above to bound the val- 
ues of A and M (or a) for different n taking /c = 1 for 
simplicity. The constraints for n = 4, 6, 10 with k = 1 
are shown in Fig. [TOl In these figures we have included 
all the bounds coming from the Eot-Wash experiment, 
as well as those coming from Casimir force searches. We 
also include the bounds (labeled Irvine) coming from an- 
other search for Yukawa forces. In general, the larger n 
(and /c), the larger the region of allowed parameter space. 
This is the case because, in a fixed density background. 



4)(2 + fc) 2 + fc 

(jn+k and 



the chameleon mass, mc, scales as M 
therefore rUc increases with n and k since the exponents 
are monotonous functions of n and k. The larger rric 
is, in a given background, the stronger the chameleon 
mechanism, and a stronger chameleon mechanism tends 
to lead to looser constraints. The chameleon mechanism 
also becomes stronger in the limits M ^ or a ^ oo, 
and all of the constraints are more easily satisfied in these 
limits. The interesting region of the parameter space is 
when M ~ M^e and \ 1. When A is very small, 
the chameleon mechanism is so weak that, in all cases, 
the chameleon behaves like a standard (non-chameleon) 
scalar field and the bounds depends solely on the value 
of |Mpi. It is clear that A ^ 1 (which implies 
1/3,0^ |Mpi 1) is very much allowed for a large class 
of chameleon theories. This is in agreement with what 
was found for the SCM in (5|. 



V. CONCLUSIONS 

We have studied a scalar-tensor theory with a field de- 
pendent coupling (assumed to be of the form of an inverse 
power-law) and a power-law self-interacting potential. 

Our main result is that this theory exhibits the 
chameleon mechanism as found in the original chameleon 
proposal [3]. Thus, the theory presented here is a 
chameleon field theory and many of the familiar prop- 
erties of the standard chameleon model carry over to 
this new setup. The thin-shell solutions we find do not 
have the geometrical interpretation as found in the SCM. 
Nevertheless we can define a thin-shell factor which is 
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FIG. 9. Constraints on chameleon theories coming from ex- 
perimental searches for the Casimir force. The shaded area 
shows the regions of parameter space that are allowed by the 
current data. The solid horizontal black lines indicate the 
cases where M and a take 'natural values'. The solid vertical 
lines show when — Hq. The dashed black line indicates 
when 1/3, |Mpi = 1 for pc — 0{lg/cm^). The amount of 
allowed parameter space increases with n. 
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FIG. 10. Combined constraints on chameleon theories. The 
shaded area shows the regions of parameter space that are 
allowed by the current data. The solid horizontal black lines 
indicate the cases where M and a take 'natural values'. The 
solid vertical lines show when — Hq. The dashed black 
line indicates when 1/3,^^ |Mpi = 1 for pc = 0{lg/cm^). The 
amount of allowed parameter space increases with n. 
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analogous to the thin-shell factor in the SCM and which 
describes the suppression of the fifth-force. 

The effective coupling in the thin-shell solutions are of 
the same form as found in the SCM when using a power- 
law potential, compare ([52|) with [5, Eq. 24]. This result 
is not surprising since the self-interactions are dominating 
when we have a thin-shell. 

If we look at the bounds computed here we see that 
the natural values M = M^e (c^ = ^r) together with 
A ~ 1 are ruled out by the Eot-Wash experiment for 
n < 10. But we have shown that there exists a large re- 
gion, M < Mde, in parameter space which is allowed by 
experiments and in which |Mpi ^ 1. These results 
are equivalent to what has been found in the SCM, and 
are due to the thin-shell effect. 

Assuming that the scalar field plays the role of dark en- 
ergy, we need to fine-tune the mass-scale in the coupling 
sector, namely, we have to demand that ~ Hq. It 
should be noted that even though this mass scale has an 
unnaturally small value, it is not strictly this scale which 
determines the coupling strength to gravity: it is given 
by |/3,0Q |Mpi which can be of order unity or larger. If 
we redefine the field by (p!9|) then this fine-tuning can be 



removed and the resulting action has only one fine-tuned 
mass-scale. 

In cosmology the field is well behaved and can act as a 
dark-energy field causing the late time acceleration of the 
universe in the same manner as the SCM. The evolution 
of the density parameters when the field is slow rolling 
along the attractor is very close to that of ACDM. The 
cosmological properties of this model are left for future 
work. 
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